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Abstract. The calculation of the conductance of ballistic rings requires a 
theory that goes well beyond the Kubo-Drude formula. Assuming "mesoscopic" 
circumstance of very weak environmental relaxation, the conductance is much 
smaller compared with the naive expectation. Namely, the electro-motive- 
force induces an energy absorption with a rate that depends crucially on the 
possibility to make connected sequences of transitions. Thus the calculation 
of the mesoscopic conductance is similar to solving a percolation problem. 
The "percolation" is in energy space rather than in real space. Non-universal 
structures and sparsity of the perturbation matrix cannot be ignored. The latter 
are implied by lack of quantum-chaos ergodicity in ring shaped ballistic devices. 



1. Introduction 

Closed mesoscopic rings are of great interest [1-13]. For such devices the relation 
between the conductance and the internal dynamics is understood much less [HI [15] 
than for open systems. First measurements of the conductance of closed mesoscopic 
rings have been reported more than a decade ago 13J . In a typical experiment 
a collection of mesoscopic rings is driven by a time dependent magnetic flux <!>(£) 
which creates an electro-motive- force (EMF) — $ in each ring. Assuming that Ohm's 
law applies, the induced current is X = — (?$ and consequently the rate of energy 
absorption is given by Joule's law as 

W = Rate of energy absorption — G$ 2 (1) 

where G is called the conductance For diffusive rings the Kubo formalism leads to 
the Drude formula for G. A major challenge in past studies was to calculate the weak 
localization corrections to the Drude result, taking into account the level statistics and 
the type of occupation [12] . It should be clear that these corrections do not challenge 
the leading order Kubo-Drude result. 

Theoretical challenge: It is just natural to ask [TU [TS] what happens to the 
Drude result if the disorder becomes weak (ballistic case) or strong (Anderson lo- 
calization case). In both cases the individual eigenfunctions become non ergodic: a 
typical eigenfunction do not fill the whole accessible phase space. In the ballistic case 
a typical eigenfunction is not ergodic over the open modes in momentum space, while 

\ Present address: Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 
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§ The terminology of this paper, and in particular our notion of "conductance" are the same as in 
the theoretical review |12j and in the experimental work (13] . 



The mesoscopic conductance of ballistic rings 



2 



in the strong localization case it is not ergodic over the ring in real space. 

Major observation: Lack of quantum ergodicity implies that the perturbation 
matrix is very structured and/or sparse. Consequently the calculation of G requires 
a non-trivial extension of linear response theory (LRT). Such extension has been pro- 
posed in Ref . [15] and later termed "semi linear response theory" (SLRT) [T6J [17] . As 
in the standard derivation of the Kubo formula, also within the framework of SLRT, 
the leading mechanism for absorption is assumed to be Fermi-golden-rule (FGR) tran- 
sitions. These are proportional to the squared matrix elements \T nm \ 2 of the current 
operator. Still, the theory of [15] does not lead to the Kubo formula. This is because 
the rate of absorption depends crucially on the possibility to make connected sequences 
of transitions, and it is greatly reduced by the presence of bottlenecks. It is implied 
that both the structure of the |X nm | 2 band profile and its sparsity play a major role 
in the calculation of G. 

SLRT and beyond: Within SLRT it is assumed that the transitions between 
levels are given by Fermi-golden-rule, but a resistor network analogy [18] is used in 
order to calculate the overall absorption. The calculation of the energy absorption in 
Eq.JT]) is somewhat similar to solving a percolation problem. The "percolation" is in 
energy space rather than in real space. A recent study [TH] suggests a way to go be- 
yond the FGR approximation. If the results of Ref [19] could be extended beyond the 
diffusive regime [20] it would be possible to extend SLRT into the non-linear regime. 

Scope: In a follow up work [21] its is demonstrated that for a very strong disorder 
SLRT leads naturally to the resistor network "hopping" picture [H] [S3] , from which 
Mott's variable-range-hopping approximation [24, 25 can be derived. In the present 
work we apply SLRT to the other extreme case, of having a ballistic ring. It should 
be appreciated that the generalized resistor network picture of SLRT provides a firm 
unified framework for the calculation of the mesoscopic conductance: The same recipe 
is used in both extreme cases without the need to rely on an ad-hock phenomenology. 

Ballistic rings: In the case of a diffusive rings the mean free path I is assumed 
to be much less than the perimeter L of the ring. In the case if ballistic ring we have 
the opposite situation I ^> L. One way to model a ballistic ring is to consider very 
weak disorder (Fig. la). Another possibility is to consider a clean ring with a single 
scatterer (Fig. lb) or a small deformation (Fig. lc). In the latter cases it is natural to 
characterize the scattering region by its total transmission gx ~ 1. Consequently the 
mean free path is 

I « - - L (2) 
1 — 9T 

The assumption of ballistic motion implies that "quantum chaos" considerations are 
important. Consequently one should go well beyond the conventional random wave 
picture of Mott [24[[25], and beyond the standard Random Matrix Theory analysis. 
In fact we are going to explain that in this limit the eigenfunctions are non-ergodic, 
and hence the perturbation matrix \I n m\ 2 is structured and sparse. 
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Outline: The purpose of this paper is to discuss the modeling of ballistic rings; 
to clarify the procedure which is involved in the calculation of the conductance; and 
to analyze the simplest example. In sections 2 we distinguish between 

• Disordered rings (e.g. Anderson model) 

• Chaotic rings (e.g. billiard systems) 

• Network models (also known as "graphs" ) 

These are illustrated in Fig. 1. In particular we motivate the analysis of a simple 
prototype network model for a multimode ring. This model has all the essential 
ingredients to demonstrate the major theme of this paper. 

In sections 3-8 we elaborate on the procedure which is used for the calculation of 
the conductance. This is not merely a technical issue, since new concepts [15j [17] are 
involved. We make a distinction between: 

• The Landauer result for the open device 

• The classical Drude result for a closed ring 

• The (quantum) spectroscopic conductance of a ring 

• The (quantum) mesoscopic conductance of a ring 

Here "classical" as opposed to "quantum" should be understood in the sense of 
Boltzmann picture. In the "classical" case the interference within the arm of the 
ring is ignored, while both the Fermi statistics and the single scattering events are 
treated properly. For the prototype model we get the following simple results: Given 
that the device has M. open modes and its total transmission is qt < 1, the Landauer 
conductance is 

G Landaucr = ^-M g T (3) 
The Drude result for the closed ring is 

G Drude = ^-M (4) 
Awn 1 — gT 

and the associated quantum results are 

9t 



e 2 

G spBC ~ — -M x minimum 
27m 



1 - gr 



,M 



(5) 



p 2 



Gmcso = TTT 2M g»~ ( 6 ) 

27m 

The calculation of g mcso involves a complicated coarse graining procedure that we 
discuss in section 8. The spectroscopic result G apoc describes via Eq. (JTJ) both the 
initial (transient) rate and also the long time (steady state) rate of energy absorption, 
provided the environment provides a strong relaxation mechanism. The mesoscopic 
result describes the (slower) long time rate of energy absorption if the environmentally 
induced relaxation is weak. See [15] for an extended quantitative discussion. 

The outcome for G mcso may differ by orders of magnitude from the conventional 
Kubo-Drude result. The calculation procedure implies that G mcso < G spoc < G Drudc . In 
the last part of this paper (sections 9-12) we demonstrate this point via the analysis of 
the prototype model. Our numerical results, whose preliminary version were reported 
in [37j, suggest that typically G mcso < G Landauor . The results of the calculation are 
contrasted with those of the conventional Kubo approach, and their robustness is 
discussed. 



The mesoscopic conductance of ballistic rings 



4 



2. Modeling 

A simple model for a ballistic ring can be either of the "disordered type" or of the 
"chaotic type" . Let us visualize the disordered potential as arising from a set of 
scatterers which are distributed all over the ring (Fig la). Depending on the scattering 
cross section of the individual scatterers we can have £ -C L for strong disorder or 
I ^S> L for weak disorder, where I is the mean free path for velocity randomization, 
and L is the perimeter of the ring. 

Another possibility is not to make all the scatterers smaller, but rather to dilute 
them. Eventually we may have a chaotic ring where the scattering is induced by a 
single scatterer (Fig lb). For example the scatterer can be a disc or a semi-disc as in 
Fig 2. These variations of Sinai billiard (billiard with convex wall elements) are known 
to be chaotic. It is important to remember that "chaos" means that complicated 
ergodic classical dynamic is generated by a simple Hamiltonian {no disorder!). 

In our view chaotic rings are more interesting for various reasons. Ballistic devices 
are state-of-art in mesoscopic experiments. For example it is quite common to fabricate 
Aharonov-Bohm devices. In such devices it is possible to induce local deformation of 
the potential by means of a gate voltage. Hence one has a full control over the 
amount of scattering. Also from theoretical point of view it is nice to have a well 
defined scattering region: This allows to use the powerful S-matrix point of view that 
has been initiated by Landauer. In particular we can ask what is the conductance of 
a device depending on whether it is integrated in an open geometry as in Fig. le or 
in a closed geometry as in Fig. Id. We believe that "chaos" and "disorder" lead to 
similar physics in the present context, but this claim goes beyond the scope of the 
present paper. 

A multimode ring can be visualized as a waveguide of length L and width W . 
In such case the number of open modes is Ai cx (k F W) 1 where d — 2,3 is the 
dimensionality. We label the modes as 



The scattering arise due to some bump or some deformation of the boundary, and can 
be described by an 2A4 x 2A4 scattering matrix S. For the semi-disc model analytical 
complicated expressions are available [29, 30J . The "classical" transitions probability 
matrix g is obtained by squaring the absolute values of the S matrix elements. 

Disregarding the closed channels, the ballistic ring is described as a set of M. open 
modes, and a small scattering region that is characterized by its total transmission gx- 
Optionally the ballistic ring can be regarded as a network: Each bond corresponds to 
an open mode. Let us consider the simplest model where the scattering is the same 
for an incident particle that comes from the left or from the right: 



where g R is the reflection matrix and g T is the transmission matrix. The simplest 
model that one can imagine is with 



a = mode index = 1,2,...,^ 



(7) 




(8) 




(9) 
(10) 



[g T ]a,b = (1 - Me 2 )S a , b 



such that the total transmission is 



g T = 1 - Me 2 



(11) 
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Such "classical" transitions probability matrix can arise if we take the S matrix as 
_ /e expj^r^ ) Vl-MeH a , b \ 

{Vl-Me 2 S a , b - e exp(-i2ir$)) [LZ) 
There a lot of simplifications that were involved in construction this S matrix. 

• The forward scattering is to the same mode only 

• The back scattering is "isotropic" 

• The scattering is energy independent 

• The scattering phases are not random 

One can wonder whether this S matrix still qualifies as 'generic', or maybe the model 
is over-simplistic. In order to illuminate this point let us look at the Sinai billiard 
models of Fig. 2. These models are fully qualified as "quantum chaos" systems. One 
observes that the specific g matrix of Eq.([9lfT0| is inspired by that of Fig. 2a. In 
this billiard an incident particle is equally likely to be scattered to any mode in the 
backward direction, but the forward scattering is only to the same mode (same angle) . 
As for the phases: we already have explained that our interest is not in disordered ring, 
but rather in chaotic one. Therefore to have random phases in the S matrix is not an 
essential feature of the model. The phases are effectively randomized simply because 
the wavenumber k a is different in each mode. Optionally, it is more convenient to 
assume that all the k a are equal, and instead to have bonds of different lengths L a . 
This provides the required phase randomization. 

Generality: Though the arguments above strongly suggest that the simplified 
network model is generic, we were careful to verify [28 that indeed all the results 
that we find are also applicable in the case of the Sinai-type system of Fig. 2b where 
the scatterer is a semi-disc. We further discuss the robustness of the results in the 
concluding section. Sinai-type billiards are recognized as generic chaotic systems. The 
reason for preferring the network model as the leading example in the present paper 
is both pedagogical and practical: The mathematics is much simpler, and the quality 
of the numerics is much better. 



3. The classical Kubo formula and Drude 



The Fluctuation-Dissipation version of the Kubo formula expresses the conductance G 
as an integral over the current-current correlation function: 

G = £ F x i f (l(t)l(0))dt (13) 

The density of states at the Fermi energy is 

g F = GeometricFactor x M — (14) 

nnv F 

where v F is the Fermi velocity. The geometricfactor depends on the dimensionality 
d = 1, 2, 3. The current observable is defined as follows: 

T = evS(x — xq) (15) 

where v and x are the velocity and the position observables respectively. In the 
quantum case a symmetrization of this expression is required, so as to obtain an 
hermitian operator. The section through which the current is measured is arbitrary 
and we simply take xq = +0, namely just to the right of the scattering region. 
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Optionally, if there were not the "black box" region of the S matrix, one could average 
over Xq, leading to X = (e/L)v. 

The term "classical" Kubo formula implies in this context that the current-current 
correlation function is evaluated classically, ignoring quantum interference. In the 
case of hard chaos system this correlation function decays exponentially. The Drude 
expression is the simplest classical approximation: 

2 



(T(t)I(O)) = i 



exp 



Substitution into the Kubo formula leads to 

e 2 £ 
2nn L 



(16) 



(17) 



where we have dropped the d dependent prefactor which equals 1 for networks [d = 1). 
In the case of a ballistic ring with a restricted scattering region (as in Fig. Id) it is 
more convenient to characterize the device by its total transmission < < 1 instead 
of the mean free path. The envelope of the current-current correlation function is 



\2g T -l\ 



<=> exp 



1*1 



(18) 



With the identification of grounds as t/(L/v F ) we deduce that for qt ~ 1 the mean 
free path is £ « L/(l — cjt)- A more detailed analysis 14] leads to the result 

c 2 



G Drud 



2irh 

a,b 



[2<7 T /U 



(19) 



If the device were opened as in Fig. le we could ignore the multiple rounds. In such 
case we would obtain 

G Landaucr = St Y\a T ] a ,b (20) 



— Y\9 T ] 

a.b 



Both results, the Drude result for the closed device and the Landauer result for the 
corresponding open geometry, are "classical" in the sense of Boltzmann, which means 
that they depend only on g. 

Let us see what do we get for G Landauor and for G Drudc in the case of the prototype 
system that we have defined in section 2. The calculation of G Land auc r is trivial and 
leads to Eq. |3|) . The calculation of G Drudo is more complicated since it involves matrix 
inversion. Still g is sufficiently simple to allow a straightforward calculation that leads 
to Eq.([4]). The rest of this section is devoted to the details of this calculation. 

We write gr = t 2 1 and gn = e 2 Y, where r and e are defined via 
2 = 1 — A4e 2 , and where we have introduced the following MxM matrices 



9t 



\ 



and T 



\ 



V- / v- 

Note that the two matrices commute. Using these notations we get 
2g T 2r 2 1 



(21) 



1 - 9t + 9r 1 - t 2 l 



2t 2 
1 - r 2 1 



1 



2r 2 



cc L 



(1-r 2 ) 2 



(22) 
(23) 
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where we have defined the normalized column vector 

c Q = -=L= a = l,2,---M (24) 

V 1 — T 

and we have used the identity 

1 1,1, 

= 1 — cc* = 1 -cc* (25 

l + cc f l + c*c 2 y ' 

Observing that J2 a b ~^- ab ~ M. and J2 a b "^ ab = -M- 2 we § et tne desired result Eq.(|4]). 



4. The quantum Kubo formula and beyond 

Our objective is to find the conductance of the closed ring in circumstances such that 
the motion inside the ring is coherent (quantum interference within the bonds is not 
ignored) . The calculation is done using the quantum version of Eq. (|13p which involves 
the matrix elements T nm of the current operator: 

G = nh gl x ((\I nm \ 2 )) (26) 

This equation would be the traditional Kubo formula if ((...)) stood for a simple 
algebraic average over near diagonal matrix elements at the energy range of interest. 
By near diagonal elements we mean \E n — E m \ < T, where T is level broadening 
parameter. The levels of the system are effectively "broadened" due to the non- 
adiabaticity of the driving [31j or due to the interaction with the noisy environment [5] . 
In what follows we assume 

A < r < A fc (27) 

where A = l/g F is the mean level spacing, and A;, = Tthv F /L is the Thouless energy. 
(Note that A;,/A = A4). Contrary to the naive expectation it has been argued in [T5] 
that depending on the physical circumstances the definition of ((...)) may involve a 
more complicated coarse graining procedure. Consequently the result for G may differ 
by orders of magnitude from the traditional Kubo-Drude result. We shall discuss this 
key observation in later sections. 

For a network system g F = ML/(irhv F ). Furthermore it is convenient to define a 
scaled matrix I nm via the relation 

Inm = -i(ev F /L)I nm (28) 

so as to deal with real dimensionless quantities. Thus we re- write Eq. (f2"6")) as: 

G = 7p-r 2M 2 g (29) 
27m 

where g = ((\I nm \ 2 )). In later sections we shall discuss the recipe for the g calculation. 
It is important to realize (see next section) that g < 1. This implies a quantum 
mechanical bound on G. 



5. The quantum bound on G 

We write the channel wavefunctions as ^ a (x) — A a sm(kx + <p), where a labels the 
modes. In our simplified network model the modes are re-interpreted as bonds, and 
we assume that the wavevector k = (2mE) 1 ^ 2 is the same for all bonds. For the matrix 
elements of X we have the expression 

I nm » -*ev F ]T Bin(¥#° " vi m) ) (30) 

a 
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were the approximation takes into account that our interest is in the couplings between 
levels with k n ~ k rn m k F . From this expression we deduce that the scaled matrix 
elements of Eq. ES]) are bounded as follows: 



Inm < E^A^A^ < 1 ( 31 ) 

a 

Therefore, irrespective of the details of the averaging or coarse graining procedure, it 
is clear that g < 1 as stated at the end of the previous section. Consequently 

G = ^-M 2 (32) 

maximal 

In the last expression we have omitted a factor of 2. This is not a typo. We shall 
explain this point in section 12. 



6. The ergodic result for G 



The simplest hypothesis is that all the wavefunctions are ergodic random waves. This 
is in the spirit of Mott's derivation [24l [25], where it has been demonstrated that 
a random wave assumption recovers (via Eq. (|2"6f ) the Drude result. If indeed the 
wavefunctions were spread equally over all the bonds, it would imply |A a | 2 ~ 2/ (ML). 
If this were true we would get 

2 

1 



Unr. 



2M 



This would imply that the conductance of a ballistic (chaotic) ring is 

„2 



G 



-M 



(33) 



(34) 



xgodic 2ttH 

We would like to argue that this result is wrong. Moreover, it must be wrong. The 
result is wrong because the eigenfunction of a ballistic ring are not ergodic. This 
we discuss in section 10. Furthermore, the result must be wrong because it violates 
quantum-classical correspondence, which we discuss in the next section. 



7. The quantum conductance and Drude 

In this section we define the distinction between mesoscopic and spectroscopic 
conductance and further discuss the latter. In the next section we elaborate 
on the calculation procedure of both. We would like to clarify in advance that 
the spectroscopic conductance is the outcome of the traditional Kubo calculation. 
Moreover, it is only the spectroscopic conductance which obeys quantum-classical 
correspondence considerations. 

Mesoscopic conductance: If the environmentally induced relaxation can be 
neglected, the rate of energy absorption depends on having connected sequences of 
transitions between levels j!5j . In the next section we explain the proper procedure 
for calculating the conductance in such circumstances. The result that comes out from 
such calculation is what we call the mesoscopic conductance G mosa . 

Spectroscopic conductance: Within the framework of linear response theory, 
it is assumed that the EMF-induced transitions are very slow compared with the 
environmentally induced relaxation. Then one can argue that Eq. (|26p is valid with 
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((|^nm| 2 )) interpreted as an algebraic average over the matrix elements. This is what 
we called in section 4 the traditional Kubo formula. Optionally, if applicable, one 
may perform an algebraic average over realizations of disorder. The latter is a very 
common procedure in diagrammatic calculations. The outcome of the (traditional) 
calculation is what we call the spectroscopic conductance G apcc . For further discussion 
of the conditions that justify a "spectroscopic" calculation see Ref. [15]. 

The spectroscopic conductance is not very sensitive to T. In fact the T dependence 
of the result is nothing else but the weak localization correction |12) . It scales like 
A/F for diffusive rings, where A is the mean level spacing (note Eq . ([27]) ) . 

Disregarding weak localization corrections it can be argued [15] that the 
obtained result for G 5poc is G Drudo provided some reasonable quantum-to-classical 
correspondence conditions (see below) are satisfied. It follows that the ergodic 
hypothesis of the previous section cannot be correct, because G Drudo is definitely not 
bounded by the number of open modes - it can be much larger. 

The necessary condition for quantum-classical correspondence can be deducted 
by taking into account the quantum bound of section 5. As gx becomes closer to 1, the 
Drude expression diverges. Quantum-to-classical correspondence is feasible provided 
the quantum bound is not exceeded: 

<C M (35) 

1 - 9t 

This can be re-phrased as 

- < M (36) 

or as 

td < t H (37) 

where t c i — £/v F is the ballistic time, and tn = M. x (L/v F ) is the Heisenberg time (the 
time to resolve the quantized energy levels). In order to establish quantum-to-classical 
correspondence in a constructive manner one should express the Kubo formula using 
Green functions, leading to a double summation over paths. Then one should argue 
that energy averaging justify the use of the diagonal approximations. The procedure 
is the same as in [32]. Needless to say that Eg. (1341 is not consistent with this argued 
correspondence and therefore a-priori must be wrong. 

8. The calculation of G 

The SLRT recipe for the calculation of G is implied by the following statements: 

[A] The transitions rates between levels {w nm ) are given by the FGR. 

[B] The diffusion in energy (D) is given by a resistor network calculation. 

[C] The diffusion-dissipation relation at low temperatures is W — g F D. 

[D] The expression for the conductance G is identified via Eq.([T]). 

In this section we give all the relevant details for using this recipe. The Hamiltonian of 
the ring in the adiabatic basis is TL i— ► E n S nm +W nm where W nm = i&KL nm / (E n — E m ) , 
and — <i> is the EMF. The FGR transition rate between level n and level m is 

2ir 

w nm = —8{E n ~E m )\W nm \ 2 (38) 
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Since we are dealing with a closed system one should take explicitly into account the 
broadening of the delta function: 

E„. — E„ 



o~(E n —E. n 



r 



(39) 



The normalized kernel F() reflects either the power spectrum or the non-adiabaticity of 
the driving. For the purpose of numerical demonstration we assume F(r) = exp(— 2\r\) 
as in [17] , The level broadening T is identical with F of Ref.[15l [31] and with fruiQ of 
Ref.|17j. As in the conventional derivation of linear response theory [T^], also here 
we regard T as a free parameter in the theory. It is convenient to use dimensionless 
quantities, so we re- write Eq. (j38|) as: 

„2 



irn 



where the dimensionless transition rates are 



2g 



-3 



If 

F 



E n — E„ 



(40) 



(41) 



(E n - E m f 

In practice we could make the approximation (E n — E m )/A sa (n— m), that 
underestimates exceptionally large couplings between almost degenerated levels. Such 
an approximation would not be reflected in the G mcso calculation (see below), because 
the latter is determined by the bottlenecks. 

The FGR transitions between levels lead to diffusion in energy space. We would 
like to calculate the coarse grained diffusion coefficient D without assuming that all 
the |/n m | 2 are comparable. For this purpose it is useful to exploit the following resistor 
network analogy |17j : 

w- 1 

d- 1 

In dimensionless units • 

units D is denoted as g and it is defined via the following equation: 

„2 



• resistor between node n and node m 

resistivity of the network 
is denoted as g„ m as defined via Eq. (|40 



D = g- 



Tit) 



M z g<$> 2 



(42) 
(43) 

In dimensionless 
(44) 



The extra q~ 2 factor compared with Eq. (|4"0")) arise because the resistivity D^ 1 is 
calculated per unit "energy length" while the scaled resistivity g" 1 is per unit site. 

A standard numerical procedure is used for extracting g for a given resistor 
network g nm . The steps are as follows: (i) Cut an N site segment out of the network 
(Fig. 3). (ii) Define a vector J n (n = 1..N) whose elements are all zero except the 
first and the last that equal J\ = +J and Jjy = — J- (iii) Solve the matrix equation 

Jn = J2 gnm ( F ™ ~ Vm ) ( 45 ) 
m 

This equation should be solved for V n . In practice it is easier to write this equation 
as J = -gV where g nm = g nm - 6 nrn J2 m > Znm' ■ The difference V\ - V N is obviously 
proportional to the injected J. (iv) Find the overall resistance of the truncated 
network g^ = J/(Vn — Vi). And finally: (v) Define the resistivity as g _1 = gJ^/N. 
For a locally homogeneous network it has been argued in Ref.[15] that g can be 
obtained via an harmonic average: 

-l 



1 

N 



N 

E 



(46) 
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The internal sum reflects addition of resistors in parallel, while the harmonic average 
reflects addition of resistors in series. This should be contrasted with the algebraic 
average which is used in order to calculate the spectroscopic result: 



1 N 

-Y 



(47) 



It is a simple exercise to verify that if all the matrix elements are the same, say 
|-fnm| 2 = c 2 , then g moso = g 3pcc = a 2 . But if the matrix is structured or sparse then 
gmeso is much smaller compared with g apcc . Schematically we write in both cases 

g = «l^™| 2 » (48) 
It should be clear that in both cases (spectroscopic, mesoscopic) the "averaging" 
requires the specification of the smoothing scale T as implied by Eq. (|4"Tj) . It is also 
clear that the mesoscopic result is much more sensitive to the value of T. Unlike 
the spectroscopic result where the dependence on T is merely a "weak localization 
correction" [12] , in the case of the mesoscopic result the dependence on T is a leading 
order effect. 

The diffusion-dissipation relation states that the rate of energy absorption is 
W = g F D. Then it is implied by Eq.([T]) that the conductance G is given by Eg. (|2"9"j) . 
The procedure above can be summarized by saying that g can be calculated from |/„ m | 2 
via an appropriate "averaging procedure" . The appropriate averaging procedure is 
algebraic (Eq. pT)) ) in the case of the spectroscopic conductance. The appropriate 
averaging procedure is harmonic-type (as discussed above) in the case of mesoscopic 
conductance. 

9. The eigenstates of the network model 

The network model that we have presented in section 2 is defined in terms of the 
scattering matrix S u , and the free propagation matrix S Wl 
(e e X p(i2n% ) VT^M?6 a , b \ 
* D \y/l -Me*6 a , b -eexp(-i2Tr%)J 

q ^ikL 



Sw ~ { e* kL ° 6 ab a ) (50) 

The wavefunction can be written as 

M 

\i>) .— > Y,{ A ^ eik{x ~ La)+A ^ e ~ ikx ) ®i a >- ( 51 ) 

a=l 

The set of amplitudes Al and An that can be arranged as a column vector of 
length 2M.. The linear equation for the eigenstates is 

£)-*.«.(£) m 

and the associated secular equation for the eigenvalues is 

det[S w S D -l] = (53) 

In the absence of driving we have time reversal symmetry, and the unperturbed 
eigenfunctions can be chosen as real (see appendix A): 

M 

\ip) i — ► J2A a sm(kx + tp a ) <g> \a). (54) 
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The wavefunction is normalized as 



E 

a=l 



A 2 a sm 2 (kx + ipa) dx = 1 (55) 



o 



which implies 



M 

Ey^-i (56) 

a=l 

For a given gr we can find numerically the eigenvalues and the eigenstates, thus 
obtaining a table 

(k„, <p£\AW) n = level index (57) 
For the numerical study we have chosen a network system consisting of M. — 50 bonds. 
The length of each bond is randomly selected in the range L a = 1 ± 0.1. We select 
the eigenvalues with k n <~ 2000. The numerical results over the whole range of gr 
values are presented in Figs. 4-7. In the following sections we discuss and analyze 
these results. 

It is of course possible to determine analytically what are the eigenvalues and the 
eigenstates in the — > 1 limit. The combined scattering matrix is 

^w^-^^+i^) T e ikL *5 a , b ) (b *> 

where we use the notation r = (1 — Me 2 ) 1 / 2 . For gr = 1 this matrix becomes diagonal. 
Then it has M. distinct eigenvalues, each doubly degenerate. We are interested in the 
non-degenerate case in the limit e — > 0. The eigenstates are still localized each in 
a single a bond, but the degeneracy is lifted. Within the framework of degenerate 
perturbation theory we have to diagonalize the 2x2 matrix 



(59) 



whose eigenvalue are determined by the associated secular equation 

(e 2 + T 2 )c 2lkL « ~ 2re lkL « + 1 = (60) 
Hence we get the following approximations 

2ir X integer ± L_ e J (61) 

We have verified that the numerical results of Fig. 4 and Fig. 7 agree with these 
estimates. We note that for hard wall scatterer each <p a would become either or tt/2 
in the gx — * 1 limit. 

10. The non-ergodicity of the eigenfunctions 

In Fig. 5 we display images of the column vectors for two representative values 
of gx so as to illustrate the crossover from localized to ergodic wavefunctions. Each 
eigen-function can be characterized by its participation ratio: 



PR 



(63) 
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This constitutes a measure for the ergodicity of the eigen-functions. By this definition 

PR ~ { ^ a s ^ n S^ e bond localized state ,g^, 

~ \ M. for a uniformly distributed state 

We distinguish between 3 regimes depending on the value of the total transmission <?r, 

• The trivial ballistic regime (1 — gx) *C 1/M for which PR ~ 1 

• The non-trivial ballistic regime \j M. <C (1 — <?t) <C 1. 

• The non-ballistic regime where gx is not close to 1 and PR ~ M. 

In the trivial ballistic regime the eigenstates are like those of a reflection-less ring with 
uncoupled modes, hence PR ~ 1. Once (1 — g T ) becomes larger compared with 1/M 
first order perturbation theory breaks down, and the mixing of the levels is described 
by a Wigner Lorentzian. The analysis is completely analogous to that of the single 
mode case of Ref. [15] . leading to PR oc (1 — gr) x M.- For gx values that are not close 
to 1 the eigen-functions become ergodic with PR ~ M.. From RMT we expect [33] 
PR ~ A satisfactory global fit, that works well within the non-trivial ballistic 

regime is (Fig. 8): 

PR w 1 + ^(1 — gr)M (65) 
o 

Our interest is focused in the non-trivial ballistic regime 1/A4 <C (1 — <?t) *C 1, where 
we have strong mixing of levels (PR ^ 1), but still the mean free path I L/(l — gx) 
is very large compared with the ring's perimeter (t L). In view of the discussion 
in section 6, it is important to realize that in this regime we do not have "quantum 
chaos" ergodicity. Rather we have PR <SC M meaning that the wavefunctions occupy 
only a small fraction of the classically accessible phase space. 

11. The calculation of matrix elements 

Given a set of eigenstates, it is straightforward to calculate the matrix elements of the 
current operator (Figs. 9-12). We recall that the scaled matrix elements are 

Inm « E ^4 n) 4 m) sin(^ - ^ m) ) (66) 

a 

with the associated upper bound 

a 

For n = m we have I nm = 1 due to normalization, and I nm = due to time reversal 
symmetry. From now on we are interested in n ^ m. There are several extreme cases 
that allow simple estimates: 

!0 for pair of states localized on different bonds 
1 for pair of states localized on the same bond 
1 for pair of ergodic states 

If we take the phases into account we get 

for pair of states localized on different bonds 

1 for pair of nearly degenerated states on the same bond 
1/(2M) for pair of uncorrelated ergodic states 
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We have already explained in section 6 that the "ergodic" hypothesis is wrong in 
the ballistic case. It should be clear that the small PR of the eigenfunctions implies 
sparsity of I nm : The matrix elements are very small for any pair of states that are 
localized on different sets of bonds. This observation is demonstrated in Figs. 9-12. 
As the reflection 1 — gx is increased, more and more elements become non- negligible, 
and the matrix becomes less structured and less sparse. 

12. Numerical results for the conductance 

Once we have the matrix elements |/ ram | 2 we can calculate G apcc using the algebraic 
average recipe Eq. ([47|) . We can also calculate G mcso using either the resistor network 
procedure or the harmonic average approximation Eq. (|46p . Fig. 13 displays the results 
for an M. = 50 network model. The conductance goes to zero for both gx — > and 
qt — > 1. The dependence on T is plotted in Fig. 14. The rough accuracy of the 
harmonic average has been verified (not displayed). 

The dependence of G on the smoothing parameter Y is easily understood if we 
keep in our mind the band profile which is illustrated in Fig. 10. In order to improve 
our intuition we show in Fig. 12 the average value of |/ n ,n+r| 2 for r — 1,2, 3, 4, 5 as a 
function of 1 — gr for r = 1,2, 3, 4, 5. 

It should be clear that the large r = 1 elements originate from the pairs of 
almost degenerate states that were discussed in section 9. Their contribution to the 
spectroscopic conductance is dominant. The upper bound Eq. (l32|) on G is implied 
by the upper bound on |/ n Tl+1 | 2 . It was already pointed out in section 11 that the 
maximal value \I nm \ = 1 is attained for the nearly degenerate states. The algebraic 
average with the interlacing vanishingly small couplings leads to the factor of 1/2 that 
was mentioned after Eg. (|3"2")) . To avoid miss- understanding we emphasize that this 
prefactor is model specific. 

On the other hand, the large r = 1 couplings almost do not affect the mesoscopic 
conductance. This is because they do not form connected sequences. Moreover, 
as implied by our calculation recipe, large value of T cannot help to overcome the 
bottlenecks. In order to get a classical result the environment should induce not only 
level broadening (which is like the 1/T 2 rate of pure dephasing in NMR studies), but 
also a relaxation effect (analogous to the 1/T\ rate in NMR). 
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13. Discussion 

In this paper we have studied the mesoscopic conductance of a ballistic ring with 
mean free path £ 3> L. The specific calculation has been done for a network model, 
but all its main ingredients are completely generic. Ballistic rings with I 3> L are 
not typical "quantum chaos" systems. Their eigenfunctions are not ergodic over 
the whole accessible phase space, and cannot be regarded as an extended "random 
wave" . Consequently the perturbation matrix T nm is highly structured and sparse, 
and we have to go beyond the Kubo formalism in order to calculate the mesoscopic 
conductance. 

Results vs expectations: The "averaging" over the matrix elements of the 
current operator should be done according to the appropriate prescription: algebraic 
scheme for the spectroscopic conductance G 5poc , and resistor-network scheme for the 
mesoscopic conductance G moso . The calculation procedure implies that 

G moBO < G apoc < G Drud( , (68) 

Our original naive belief, before we started with the numerical work, was that it is 
feasible to get quite large G moso , possibly of the order of G spoc . To our surprise the 
numerics has revealed that typically 

G moso < G Landauor (69) 

We have pushed our numerical verification of this statement up to M. = 450 (Fig. 14b). 
For an optimal value of T, such that G mcso is maximal, we still have G moao < G Landaucr . 
The numerical prefactor in the latter inequality appears to be roughly 3, but obviously 
we cannot establish that there is a strict limitation. Still, as far as order of magnitude 
estimates are concerned, our conjecture is that this statement is true in general. We 
did not find a mathematical argument to establish this conjecture, except the very 
simple case of a single mode ballistic ringfTHj where the calculations of G can be done 
analytically. 

Robustness: Our results are not sensitive to the details of the model. 
Disregarding the details, the eigenfunctions are doomed to be non-ergodic in mode 
space if 5t ~ 1- This by itself implies that I nm is sparse and possibly structured. 
Consequently the resistor network picture implies that G moso is much smaller compared 
with the naive expectation. We have further tested the generality of our quantitative 
statements by analyzing [28] another, more realistic model, where the ring is modeled 
as a waveguide with a semi-disc scatterer (Fig. 2b). The S matrix for this model 
is known [29[ 130] . In particular we have verified that the participation ratio of the 
eigenstates has roughly the expected dependence on gx- Indeed for both the semi- 
disc model, and our simplified network model, the participation ratio does not exhibit 
anomalous saturation as typical, say, for a "star graph" 34 . Eventually we have 
verified that the results for the mesoscopic conductance of the semi-disc model are 
similar to those that were obtained for the network model. 

Challenges: It is still an open challenge to derive an estimate for the mesoscopic 
conductance in terms of gx- It was possible to derive such an expression in the 
single mode case. There we have found that G mcso oc (1 — gr) 2 gr- In the general case 
(A4 > 1) the calculation is more complicated. We suspect that our expression for the 
participation ratio Eq. (|65|> constitutes an important step towards this goal. In any 
case we were not able to derive a reliable closed analytical expression. 

Limitations: It should be emphasized that if there is either a very effective 
relaxation or decoherence process, then the semi-linear response theory that we have 
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discussed do not apply. In the presence of strong environmental influence one can 
justify, depending on the circumstances [13], either the use of the traditional Kubo- 
Drude result, or the use of the Landauer result. 



Appendix A. Implications of time reversal symmetry 

We can decompose the eigenstate equation as follows: 

£i)= s °(t) 

£)-*•(&) ,A - 2) 

Above B^a and B Ra are the amplitudes of the outgoing waves from x = 0, while 
A^a and A Ra are the amplitudes of the ingoing waves. Conventional time-reversal- 
symmetry implies that both ip(x) and its complex-conjugate ip{ x )* satisfy the same 
Schrodingcr equation. Complex conjugation turns out the incoming wave into outgoing 
one and vice versa, and therefore 

tH(t) (A - 3) 

B* \ ( A* \ 

TD* a ) — S-w [ A* a ) (A-4) 

°Ra J \ A Ra J 

It is not difficult to see that the two sets of equations are equivalent provided 

^transposed = ^ (A 5) 

S* w p = S* w (A.6) 

If we have this (conventional) time reversal symmetry, the unperturbed eigenfunctions 
can be chosen as real in position representation: 

M 

\$)=S^A a wa.(kx + ip a )®\a). (A.7) 



where 



A a = 2\A La \ =2\A Ra \ (A.8) 
fa = - (tt + arg(A La /A Ra ) - kL a ) (A.9) 
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Fig.l: (a) A ring with disorder. The mean free path can be either I -C L for diffusive ring or I 3> L 
for ballistic ring, where L is the length of the ring, (b) A chaotic ballistic ring. Here we have a 
single scatterer. The annular region supports M open modes, (c) Another version of a chaotic ring. 
Here the scattering is due to a deformation of the boundary, (d) A chaotic ring can be regarded 
as a network. Namely, each bond corresponds to an open mode. In the numerics the lengths of the 
bonds (0.9 < L a < 1.1) are chosen in random. The scattering is described by an S matrix, (e) The 
associated open (leads) geometry which is used in order to define the S matrix and the Landauer 
conductance. 





Fig. 2: (a) upper panel: a waveguide with convex scatterer. This geometry has inspired our simple 
network model, (b) lower panel: the semi disc model. For this geometry we have some preliminary 
numerical results that will be published elsewhere. 
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Fig. 3: Within the framework of the Fermi golden rule picture the flow of the probability current in 
a multi level system is analogous to the flow of current via a resistor network. Thus the inverse of 
the course grained diffusion coefficient can be re-interpreted as the resistivity of the network. On the 
right we display a truncated segment, where + J is the current injected from one end of the network, 
while — J is the same current extracted from the other end. The injected current to all other nodes 
is zero. The resistance of each "resistors" in the network corresponds to g„m- 




Fig. 4: The eigenvalues k n within a small energy window around k ~ 2000 are shown as a function 
of the reflection 1 — g^p. We consider here a network model with M = 50 bonds. The length of each 
bond was chosen in random within 0.9 < L a < 1.1. 




10 20 30 40 50 10 20 30 40 50 

Level index (n) Level index (n) 



Fig. 5: Each culumn is a grey-level image of one eigenvector [Ai | 2 , where a = 1..M is the bond 
index. We display the eigenvectors in the range 2000 < k < 2031. Left panel: gx = 0.999. Right panel: 
9T = 0.5. 
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Fig. 6: The amplitudes |Aq | 2 with a = 1...M of one representative state (fc n ~ 2011) as a function 
of the reflection. The wavefunction is localized on a single bond for small reflection, and becomes 
ergodic for large reflection. 
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Fig. 7: The phases tp a for the same eigenstate of Fig. 6. The solid lines are the values which are 
implied by Eg. (I62I I . The crosses indicate the phases within the bond a where most of the wavefunction 
is localized. Indeed in the limit gx — > 1 this phase coincides with one of the predicted values. 
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Fig. 8: For each value of Qt we calculate the participation ratio (PR) for all the eigenstates. We display 
(as symbols) the minimum value, the maximum value, and a set of randomly chosen representative 
values. The solid line is the average PR, while the dotted line is Eq.l l65l l. The left panel is log-linear 
while the right is log-log. 
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Fig. 9: Images of [/ nm | 2 . The main diagonal is eliminated from the image. In left upper panel we 
display a relatively large representative piece for gx = 0.9. In the other panels we display zoomed 
images for qt = 0.999, 0.9, 0.5. As the reflection 1 — qt becomes larger, more elements become non- 
negligible, and the matrix becomes less structured and less sparse. 
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Fig. 10: Images of \I nm \ 2 - The main diagonal is zero due to time reversal symmetry. In left upper 
panel we display a relatively large representative piece for gx = 0.9. In the other panels we display 
zoomed images for gx = 0.999, 0.9, 0.5. As the reflection 1 — gx becomes larger, more elements become 
non-negligible, and the matrix becomes less structured and less sparse. 
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Fig. 11: The n-averaged value of 2A4\I n ,n+r\ 2 as a function of 1 — gx for r = 1, 2, 3, 4, 5. The ergodic 
value for this quantity (2M) is indicated by the solid horizontal line. We also indicate the value M 
by a dashed horizontal line. The left panel is normal scale, while the right panel is log-log. 
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Fig. 12: The n-averaged value of 2_M|/ n ,n+r| 2 as a function of 1 — gx for r = 1,2, 3, 4, 5. The ergodic 
value for this quantity (1) is indicated by the solid horizontal line. We also indicate the maximal 
value M by a dashed horizontal line. The left panel is normal scale, while the right panel is log-log. 




Fig. 13: (a) Left panel: The mesoscopic conductance G in units of e 2 /(2nh) as a function of 1 — gx- 
The curves from bottom to top arc for 7 = T/A = 1,2,3,4,5. The total number of open modes 
is M = 50. The dotted line is GLandaucr while the dashed line is GDmdc- (b) Right panel: The 
mesoscopic conductance (lower solid line) is compared with the spectroscopic conductance (upper 
solid line). Here 7 = 3. The dotted and the dashed lines are as in the left panel. 
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Fig. 14: (a) Left panel: The mesosocopic conductance as a function of 7 = T/A. The curves 
from top to bottom are for gx = 0.8,0.7,0.5. The number of open modes is hA = 50. (b) 
Right panel: The mesosocopic conductance divided by the number of modes for gx = 0.8 and 
M = 50, 100, 150, 200, 250, 300, 350, 400, 450. 



